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One approach to the understanding of fusion in cells and model membranes involves stalk for- 
mation and expansion of the hemifusion diaphragm. We predict theoretically the initiation of 
hemifusion by stalk expansion and the dynamics of mesoscopic hemifusion diaphragm expansion 
in the light of recent experiments and theory that suggested that hemifusion is driven by intra- 
membrane tension far from the fusion zone. Our predictions include a square root scaling of 
the hemifusion zone size on time as well as an estimate of the minimal tension for initiation of 
hemifusion. While a minimal amount of pressure is evidently needed for stalk formation, it is not 
necessarily required for stalk expansion. The energy required for tension induced fusion is much 
smaller than that required for pressure driven fusion. 
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I. INTRODUCTION 

Membrane hemi fusion is a possible pathway (see 
Muller et alJ ((2002) for an alternative view) to the com - 
plete fusion of membranes (jChernomordik et all Il995|) . 
Current theories associate the initiation of hemifusion 
with the formation of a contact zone between the mem- 
branes in which the two proximal monolayers are con- 
nected by a stalk-shaped neck. The stalk then expands 
and a region is formed (region C in Fig. [Q, in which 
the two distal monolayers form a single bilayer. In gen- 
eral, the energetic cost of the splay of the lipid chains in 
the stalk, prohibits its spontaneous expansion. However, 
the presence of additional, external forces (e.g. pressure, 
surface tension gradients, electrostatic effects) can lead to 
expansion of the stalk into a 'hemifusion region' and to 
the growth of this zone. Clear evidence for the existence 
of these two distinct pre-fusion stages, stalk formation 
and he mifusion, was found fo r PEG mediated fusion of 
vesicles l)Lee and LentzllT997l) . 

A recent theoretical paper l)Safran et all l200lj) sug- 
gested that the flow of lipids from region B to region 
A can be caused by an increase of the surface tension 
in region A due to the presence (in that region only) of 
additional polymer in solution. The tension gradient be- 
tween these regions induces a flow of lipids, that leads to 
the growth of region C. 

A different scenario, where hemifusion can be an 
alternative pathway to fusion wa s found in influenza 
hemagglutinin-mediated fusion llChernomordik et all 
ll998t lLeikina and Chernomordikl l200Cl). The initial lo - 
cal stalk may evolve to a fusion pore (Muller ctail l2002ft . 
or it may expand to hemifusion. In the latter case, no 
fusion occurs. 

In this paper, we predict the dynamics of the ex- 
pansion of the initial stalk and its role in the growth 
of a mesoscopic hemifusion diaphragm. The nucleation 
of a stalk by thermal fluc tuations was recently shown 
to be thermally accessible l)Kozlovskv and Kozlovl l2002t 
iMarkin and Albanesll200'2]) . A detailed description of the 
kinetics of this nucleation event (that typically describes 
the formation of a stalk of several nanometers in extent) 
is outside the scope of our work. Instead, we focus on 



estimates of the conditions that facilitate stalk expan- 
sion into hemifusion. We discuss the implications of our 
theory on biological fusion mechanisms and on in-vitro 
experiments. In addition, we predict the growth of the 
hemifusion region (e.g. from nanometers to microns) as 
a function of time and discuss the physical parameters 
that can be used to control the time scale for hemifusion. 
This dynamic part is relevant mainly to in-vitro experi- 
ments, since biological fusion events generally remain at 
the microscopic scale of the stalk. 

If hemifusion is an intermediate state of fusion then 
it is important to contrast the time scales of hemifu- 
sion diaphragm expansion and pore for mation, in order 
to det ermine the rate limiting step. IChizmadzhev et alJ 
(2000) predicted that pore expansion is exponential in 
time, with a time scale of rj m /Sp < 1 sec, where r\ m is the 
membrane viscosity and Sp is the surface tension differ- 
ence (both are estimated below). However, if pore nucle- 
ation is slow enough significant expansion of the hemifu- 
sion diaphragm can occur before pore formation. This is 
the case considered here, where we predict that the hemi- 
fusion diaphragm expands as the square root of time. 

Our theoretical model is motivate d by and consisten t 
with the experiments described by iKuhl et"afl |l996), 
where two bilayers supported on mica surfaces were 
brought into contact in the presence of a PEG-water 
solution. Hemifusion, that eventually extended over a 
distance of 50/z was observed in a time of about 10 min- 
utes, while the time it took the initial stalk to form was 
less then 3 minutes. This suggest that, at least in this 
experiment, the rate limiting step for hemifusion is the 
expansion of the fusion zone, as opposed to stalk forma- 
tion. 

This paper presents a simp le theoretical mod el relevant 
to this experimental system (IKuhl et all 11996^ , and pre- 
dicts the time dependence of hemifusion expansion. The 
overall time scale we fin d is co mparable with the mea- 
surements of IKuhl et al.l l|1996|) while the details of the 
predicted temporal dependence have yet to be tested ex- 
perimentally. 
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FIG. 1 Illustration of the experimental geometry, adopted from lKuhl et al.| Jl 
R and Rb are the inner and outer radii of region B, respectively. 



Regions A, B and C are defined in the text. 



II. PHYSICAL MODEL 

Our theoretical model is illustrated in Fig, fflthat is a 
simpli fication of the experimental system of iKuhl et alJ 
(1996) wherein two bilayers deposited on mica cylinders 
are brought together in a solution of PEG and water. 
The lipids of the distal monolayers are physisorbed on 
the mica; this fixes their lateral density. From here on in 
this paper, the term 'lipid density' relates to the lateral 
density of the proximal monolayers (see fig. LB . 

We assume that the lipids are in local equilibrium, so 
at a particular location r, the free energy per lipid (in 
the proximal monolayers) fx(f f ), does not depend on the 
lipid microstate, but only on the lipid density a(r). This 
assumption of local equilibrium is consistent with our 
results that predict an overall time scale for hemifusion 
expansion that is much larger than the local diffusion 
time of a single lipid molecule. 

The experimental system we consider is macroscopi- 
cally cylindrically symmetric and we therefore assume 
cylindrical symmetry of all the physical quantities at 
mesoscopic length scales. This is justified because all 
flows (of water and lipids) are laminar, and there are no 
mechanisms that might induce angular fluctuations or 
instabilities. 

We distinguish between three regions, illustrated in 

Fig. □ 

• Region A - where the distance d between the bi- 
layers is typically m uch larger than the polymer 
correlation length £ ijSafran et all 12001). In this 
region, the outer lipid monolayer is in contact with 
the PEG in the solution. The free energy per 
molecule in this region is given by /2(<r(r)), and is 
different (in its functional form) from the free en- 
ergy fi(a) of the monolayer in the absence of PEG. 

• Region B - where d <C £. For these values of d, 
the PEG density near the bilayers is negligible and 
our model assumes that there is no PEG in contact 



with the bilayers in this region. The free energy 
per lipid in this zone is /x(cr(r)). In addition we 
assume that the distance between the mica surfaces 
is constant (the mica surfaces in the experiment are 
deformed and flattened under pressure), and that 
this region is ring-shaped with an outer radius Rb, 
and an inner radius R. 

• Region C - the region where the distal bilayers are 
in contact. 

The bilayers are Langmuir-Blodgett deposited in wa- 
ter, without PEG, which is added later. The energy per 
lipid when the monolayers are in contact with water is 
n(cr) and the proximal monolayers are initially Langmuir- 
Blodgett deposited with a density <To that minimizes [i. 
When PEG is added, it induces an effective attraction be- 
tween the polar heads ijSafran et alll200lf) . and changes 
the functional form of the energy as function of the lipid 
density to p,{a). 

The e ffect of lipid condensation in the presence 
of PEG llBartucci et all Il996t iMaggio and Lucvl Il978t 
iTilcok and^KsherTTh^TE ) hasa been discussed in terms 
of the the dehydrati on of the bilayer by the PEG 
i|Mishima et all Il997(l . This dehydration affects the 
lipids in region A that are in microscopic proximity to 
the PEG, but has no effect on the lipids in region B. 
In Sec. IVII we demonstrate that the osmotic pressure 
induced by the PEG is too small to induce hemifusion. 
This stands in contrast to the surface tension effects that 
are the main focus of our work. 

If equilibrium could be reached, the lipid density in re- 
gion A would tend to increase in the presence of PEG. 
However, the number of lipids in the monolayers cannot 
increase to any significant degree within the time scale 
of the experiments, since the concentration of lipids in 
the bulk solution is negligible and the number of lipids 
that can be transported from region B to region A is 
much smaller than that of region A. Thus, the lipid den- 
sity is unchanged and the energy per lipid in region A is 
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now /S(cxo) > /■'('''o)) with the derivative /2'(ao) < due 
to the induced head attraction. This condensation effect 
thus leads to a negative tension in the proximal monolay- 
ers that ideally would cause them to contract in extent. 
They cannot do this without exposing the chains of the 
inner monolayers to the water and this is energetically 
prohibitive. The outer monolayers are therefore stressed 
and one way of relieving that stress is for additional lipid 
to enter this region; this will allow the local lipid density 
to increase, while still covering the original area occupied 
by the outer monolayer. 

The PEG concentration near the outer monolayers in 
region B is given by cb — CA{d/£,) 2 , where ca is the PEG 
concentration near the outer monolayers in region A, £ 
is polymer correlation length and d is the dista nce be- 
tween the bilayers in region B ijSafran et alll20fTi|) . Since 
by the definition of region B, the bilayer spacing in that 
region is small, d <C £, we have cb <C ca and the PEG 
concentration in region B is negligible; we thus take this 
concentration to be zero. The energy per lipid in region 
B is initially given by fi(o~o), where o~q is the lipid density 
in the absence of polymer. Since the free energy per lipid, 
/x, is minimized when the density a — o~o and the tension 
in region B initially vanishes, since either expansion or 
compression of the lipids will increase their energy. The 
tension gradient between regions A (initially at negative 
tension) and B (initially at zero tension) induces a flow of 
lipids from region B to region A. Since region A is much 
larger than region B, we can treat it as a reservoir, and 
assume that even though lipid is flowing from region B 
to region A, the lipid density in region A is not changed 
from its initial value of Co- The system is a dynamical 
one and the chemical potential (equivalent in our single 
component system to the free energy per lipid, n) is not 
constant in all of space at the mesoscopic scale; this re- 
sults in lipid flow and dynamics. However, since local 
equilibrium is maintained, we must have equal chemical 
potentials at any given point in the system. In particular, 
at the boundary between regions A and B, the chemical 
potentials of the lipids must be equal: n((Jb) = ft(o~o), 
where at, = <t(Rb) is the lipid density at the edge of re- 
gion B. We note that this equality of chemical potentials 
determines the lipid density at the boundary of region B, 
o~b] the functional form of the two free energies /x(cr) and 
jS(cr) are not the same, since in region A, the lipids are 
in contact with polymer. 

The initial lipid density in region B (o~q, which is the 
density at which the lipids self-assemble in water in the 
absence of polymer) is higher than the lipid density at 
the AB boundary: o-q > ov This inequality is a conse- 
quence of the fact that the tension at the boundary is 
negative, as shown in section IIVI More intuitively, the 
negative tension in region A tends to pull in additional 
lipids from the boundary region of region B into region A 
as explained above. This lipid flow reduces the lipid den- 
sity at the boundary r = Rb from cr to ov In turn, the 
reduced lipid density at the boundary of regions A and B, 
(<Tb < (Tq) induces a flow of lipids from the rest of region 



B towards the boundary. This is because the minimum 
energy state in region B is one where a = o~o > er^; thus 
lipids from the entirety of region B flow to the boundary 
in an attempt to restore the lipid density there to values 
closer to cto- This flow, in turn, reduces the lipid density 
at the boundary between regions B and C (the hemifu- 
sion region) at r = R, and lead to a negative tension that 
tends to expand region C. 

At the boundary of regions B and C, the lipid density 
is determined by a force balance between the membrane 
negative tension (arising from the lipids flowing to the 
AB boundary), that tends to expand region C, and the 
force exerted by the boundary ring around region C that 
tends to shrink it. The main contribution to the energy 
of this ring is of the tilt of the lipid tails imposed by 
the toroidal geometry. This tilt is needed in order to 
form the three-way junction of the boundary ring cross 
section while avoiding an intra- membrane void, which 
has a much higher energetic cost l)Kozlovskv and Kozlovl 
2002). The energetic cost of the tilt can be considered 
through the relate d intra-membrane strain a nd the adja- 
cent stress tensor ijHamm and Kozlovl |2000() . 

We assume that for R ^> d the energetic cost f t for 
a cross section of the BC boundary ring is independent 
of R. Thus, the ring energy is given by E r (R) = 2nRf t . 
The force per unit length that the ring exerts on region B 
of the membrane tends to shrink region C and pull region 
B in the — f direction. This force (per unit length) is 

— (-™L)r = -*f (1) 
2ttR \ dR) R ' 1 ' 

and tends to shrink the boundary ring; that is, the ex- 
pansion of region C is energetically costly. In local equi- 
librium, this force is balanced by the surface tension p, 
which may be considered as a two dimensional lateral 
lipid pressure , in region B of the monolayer that tends 
to expand the ring: 

p{R) + - = o ■ (2) 

Negative tension in region B tends to cause this region 
to contract and thus provides a force in the f direction, 
balancing the force due to the BC boundary. 

III. MONOLAYER DYNAMICS 

In this section we derive the dynamics that govern the 
expansion of the hemifusion region and predict the flow 
of lipids within the monolayer as a function of the lipid 
density and of time. 

There are three local, dissipative forces that oppose 
any lipid motion. 

• The stress, or force per unit area due to the viscos- 
ity of the water that is moved along with the lipids 
is given by r)dv w /dz, where v w is the water veloc- 
ity and r\ = 0.01 erg s/cm 3 is the viscosity of water. 



FIG. 2 Enlargement of region B (between the dashed lines) shown in Fig. Q The lipid density in region A is the initial density 
(Jo. In our model we assume a step in the density profile, so the lipid density at Rb is ab < <jq- At R it is approximated by 
go . The normalized lipid density profile p = [a — ao)/uo as a function of the radius r for R/Rb = 0.2 is plotted using Eq. 
The free energy per lipid in region B is /i(cr) which is a function of the local lipid density. In region A, the free energy per lipid 
/i(<7o) everywhere. 



The stress is of order rjv/d, where v is the lipid ve- 
locity and d is the spacing between the bilayers in 
region B. 

• The stress, or force per unit area due to the mono- 
layer viscosity is given by r^V 2 ^, where r\ m . is the 
mono layer friction coefficient (ISeifert and Langerl 
1993). For a laminar flow we estimate this stress 
as Tjmv/RB 2 ', that is, the relevant dimension is the 
size of region B in which there is monolayer flow. 

• The stress, or force per unit area that is due to 
the friction between the monolayers is given by bv, 
where b is the friction coefficient. This stress de- 
pends only on the motion of the outer relative to 
the inner monolayer where there is no flow; there is 
therefore no dependence of the length scale related 
to the geometry of the different regions. 

The friction between a DMPC monolayer and a sup- 
porting HTS 1 monolayer at T = 45°C is b = 7 • 10 6 erg 
s/cm , while for a supported OTS 2 monolayer the fric - 
tion is b = 2.9 • 1 8 erg s/cm 4 llMerkel et all Il989|) . 
The experiments of iKuhl et alJ l)l996j) were carried out 
at 25° C. It has been observed that the diffusion coef- 
ficient of a molecule in a DMPC monolayer increases 
about three folds when T is increased from 25°C to 
45°C (lHaibel et al.Lll998tlMerkel et all Il989l IVaz et all 
Il985l) . which suggest a corresponding decrease in b. In 
this work we use an estimated value of b = 10 8 erg s/cm 4 . 
For DMPC bilayers at T = 25° C the b i layer viscosity is 
r] m ~ 3 • 1CT 7 erg s/cm 2 llMerkel et allll989l) The val- 
ues relevant to the experiments of lKTrnleT'al"] l)l99fj) are 



d = 2 • KT 7 cm and R B = 5 • 1(T 3 cm. With the esti- 
mates for the stress given above, we find that the fric- 
tional force due to relative motion of the two monolayers 
is much larger than either the lipid or water viscosity 
contributions to the stress. We thus neglect these latter 
two effects and predict the dynamics for a system where 
the only relevant dissipation is due to the relative friction 
between the monolayers. 

The lipid flow is induced by the tension gradient Vp, 
and is opposed by the frictional bv. The force balance 
equation is 



Vp-bv = 



(3) 



In Appendix lAl we derive the lipid local dynamics using 
Eq. and the continuity equation. We consider the 
dynamics only to first order in the lipid density variations 
p = (a — on) /on, whi c h is k nown from experiments to be 
small. In lKuhl et alJ l|l996l) a variation of \p\ w 0.05 was 
measured. 

To first order in p the local dynamics has the form of 
a diffusion equation 



2e = 2 V 2 

dt b 



(4) 



where a — <to 3 //'(<7o) is the harmonic 'spring constant' 
of the monolayer. For a small density variation \p\ -c 1 
the surface energy cost is Sg = ^ap 2 , and the related 
tension difference is dp — ap. We note that the surface 
energy g = 074(17) is the Gibbs free energy per unit area, 
and is different then the surface tension p, which has the 
thermodynamic role of the two dimensional pressure. 
We estimate a using the phenomenological form 



1 trichlorosilanes with hexadecyl chains 

2 trichlorosilanes with octadecyl chains 



p{a) 



err/ 



(5) 
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where 7 is the effective surface tensio n of the 
hydrocarbon- water interface l)Ben-Shaull fl995). The sec- 
ond term in Eq. accounts for the (electrostatic) effec- 
tive head-group repulsion, while the first term represents 
the effective hydrocarbon-water repulsion. We note that 
this effective repulsion is smaller than the repulsion of 
the bare hydrocarbon-wat er interface, and ha s been esti- 
mated as 7 ~ 20 erg/cm 2 ijlsraelachvilil Il99l|) . 

From Eq. [5] we obtain a — 27 ~ 50 erg/cm 2 . For 
b = 10 8 erg s/cm 4 the effective 'diffusion constant' is 
5 • 10 -7 cm 2 /sec. This quantity is larger than the ac- 
tual, microscopic diffusion constant measured for free 
liquid bilayers above the gel tran sition, that are of 
the order of lO" 8 - 10" 7 cm 2 /sec llHaibel et all Il998l 
ISonnleitner et all 119991 IVaz et all Il985h . The Einstein 
relation is not applicable in our case, since the flow (that 
happens to scale like diffusion) of the lipids from the high 
to low density regions is not due to the random motion 
of the molecules, but due to the tension gradient aVp. 
Indeed, for a characteristic molecular area a = 10~ 14 
cm 2 we find that the related energy per molecule is 
era ~ lOfceT. 



IV. BOUNDARY CONDITIONS AND GLOBAL 
DYNAMICS 

The boundary conditions for the lipid density were al- 
ready discussed in section [H] and we review them here for 
convenience. The local tension equilibrium at the bound- 
ary with region A determines the local lipid density aj, 
at Rb- In Appendix lAl we show that the tension in the 
monolayer is given by 



p = a 2 p'{a) 



(6) 



Since the tension in region A is negative, from the tension 
equality at the boundary we see that p(Rb) = c& 2 //(cf,) 
is negative. Moreover, because the function p(a) has a 
minimum at <7q it is convex in a neighborhood of cfq. If 
07, is in that neighborhood, then the condition p' {ub) < 
yields that Ob < <7q. 

The boundary ring near the hemifusion region at R 
exerts a force that opposes hemifusion expansion; this 
is because the boundary energy of the hemifusion region 
is increased as this region grows. This force is locally 
balanced by the negative tension in region B where lipids 
are flowing towards region A. As lipids pass from region B 
to A the lipid density in region B decreases; the tension in 
region B, and in particular near its boundary with region 
C, becomes more negative and pulls on region C causing 
its expansion. 

The density of lipids in region B at the boundary R is 
determined from the force balance Eq. [2J Using Eq. [B] 
we may write Eq. [3 as 



P (R) 



A 

aR 



(J) 



Before the flow begins, the initial lipid density in region 
B is <7o: which implies that p — 0. For this value of the 



lipid density there is zero tension in region B, the stalk 
does not expand and hemifusion does not develop. Due 
to the tension gradient between region B and A, lipids 
flow out of region B and a negative tension is built up. If 
at a certain time the lipid density at r = R is low enough 
so that Eq. [7| is satisfied, the stalk begins to expand. 

After the flow of lipids is initiated, lipids are removed 
from region B as they flow towards region A and the lipid 
density in region B is lower than cfq. The lipid density 
in region B cannot, however, be smaller than the value 
of erf,, because when a — Ob the free energies per lipid in 
regions A and B are equal, and the flow stops. Thus, we 
require <jq > a > Ob in all of region B if there is to be flow 
and stalk expansion that leads to hemifusion. At an early 
time after the stalk formation, while the stalk does not 
expand, the lipid density in all of region B approaches the 
equilibrium density profile a(r) = Ob- Using Eq. |7| the 
condition for the stalk to begin to expand with a finite 
amount of time is: 



pb > 



ft 
aRo 



(8) 



where pb — (crb ~ &o) / &o an d Ro is the radius of the stalk. 
In our model, we consider the process for R much larger 
than the molecular size Rq that characterizes the size of 
the stalk. The tilt energy f t is in general positive. From 
Eqs. EJand|Hl for R > R we have \p{R)\ < \pb\- Since 
we consider all quantities only to first order in pb we use 
the approximation p(R) = 0. 

In Appendix[B]we use the integral continuity equation, 
that expresses the conservation the lipid number in the 
system, in order to obtain a dynamic equation for the 
hemifusion radius R. In Appendix [U] we show that the 
time scale that governs the local dynamics is much faster 
then the rate of change of R. We use an adiabatic ap- 
proximation in order to solve the dynamics. First, we fix 
R and find the asymptotic (t — > 00) lipid density profile 



' log(r/gg) \ 



(9) 



We use this density profile to obtain the dependence of 
the hemifusion radius R on the time t to find: 



(10) 



This predicts an approximately square root dependence 
of the hemifusion region size on time (with logarithmic 
correction s). The same tempor al dependence was ob- 
tained by iKumenko et all |l999) under the assumption 
of constant lateral lipid density. However, their result is 
quantitatively different from ours since they have consid- 
ered the monolayer viscosity as the main dissipative force, 
while we have showed that it is negligible compared to 
the friction b. 

From Eq. ^| we find that the time it takes the 
hemifusion region to evolve from the initial stalk ra- 
dius R = Rq <§C Rb to a final radius of R = Rb is 
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At = -bR B 2 /2ap b . With p b = -0.05 and a/b = 5 • 10~ 7 
cm 2 /sec, we predict that the time for expansion of the 
hcmifusion zone to a scale of Rb = 50 /im i s At ~ 500 
sec. T his is consistent with the experiment of lKuhl et alJ 
i|l996|) where a time of At = 600 sec was measured. 

The time At found here can also be derived (up to a nu- 
merical factor) from a simple scaling argument, that does 
not depend on the specific details of our model. As hemi- 
fusion is initiated, the tension difference between the bulk 
(at Rb) and the hemifusion front (at R) is —ap b . When 
Rb 3> R the average tension gradient is Vp « —apb/RB- 
For a fully damped flow with a friction coefficient b the 
average lipid velocity is v — Vp/6. The hemifusion front 
(BC boundary) advances with the velocity ~ v of the 
lipids near it. The time to advance a distance of Rb 
with a velocity v is Rb/v — —bRs 2 / 'apb- 



V. INITIATION OF HEMIFUSION 

The change in the monolayer surface energy due to 
the presence of PEG in region A is Sg = ao(p(<Jo) — 
/i(ero)), where /i(cro) is the free energy per lipid in the 
absence of PEG, and p,(<Jo) is the free energy per lipid 
in the presence of PEG. Since we have defined a b by the 
condition p(<Jb) — p(&o), we can expand p around its 
minimal value a = ao, and find that to lowest order in p b 
the surface energy difference Sg and the tension difference 
Sp induced by the PEG are 



SNARE 3 proteins that promote exocytosis in nerve 
synapses are thought to induce stalk expansion through 
a conformational chan ge by which the protein pull on 
the stalk to widen it l|Scales et alll200l|) . Another pos- 
sible cause for stal k expansion i s calcium ions induced 
membrane tension J Arnold H99! . We conclude from our 
theory that the latter mechanism may be more effective 
energetically. 

In section ITVl we found that in order for expansion of 
the hemifusion region to occur, the driving force due to 
the negative tension in region B must be large enough 
to overcome the tendency of the boundary of region C 
to shrink. We thus deduced that the normalized lipid 
density at Rb must obey 



From this condition, we estimate the minimum stalk ra- 
dius Rq for which the lateral tension in the monolayer 
can induce expansion. The energy of the lipid tails tilt at 
the he mifusion front is estimated bv lMarkin and Albanesl 
( 2002) as ft — 2 ■ 10~ 6 erg/cm. For the values of a and pb 
given above, we find that the mechanism described here 
is sufficient to cause hemifusion for Rq > 8nm, which is 
of the order of the typical radi us of a thermally nucleated 
stalk ijYang and HuaneL |2002|) . Note that if p b vanishes 
(that is, no polymer is present in region A) hemifusion 
will not be initiated for any finite (reasonable) stalk ra- 
dius. 



Sg = 2 a Pb 2 5 S P = ap b . (11) 

In iKuhl et ail l|l996h a change of p b « -0.05 in lipid 
density was deduced from the measured thinning of the 
bilayer. Using the value a — 50 erg/cm 2 we estimate 
Sg « 0.06 erg/cm 2 ; Sp « 2.5 erg/cm 2 . 

Initiation of stalk expansion is relevant not only to 
events of mesoscopic fusion, but also to in-vivo fusion 
events, where a fusion pore is formed soon after stalk 
expansion. In many cases of biological interest, the fu- 
sion process is regulated by fusion proteins that pro- 
mote stalk formation and expansion. One hypothe- 
sized bio-molecular mechanism that promote expansion 
is the penetration of hydrophobic fusion protein domains 
into the membrane and its subsequent destabilization 
l|Bentz and Mittal . l2000(l . The protein domains may in- 
crease the membrane surface energy by inducing an effec- 
tive attraction of the hydrophobic head groups, similar 
to the effect of PEG ijSafran et all 12001^1 : they may also 
penetrate the membrane, increasing the intra-membrane 
tension. Our theory suggests that the former mechanism, 
which work to increase in Sg, may be more effective en- 
ergetically than the latter, which increases Sp. That is, 
for a given change in lipid density, pb, a smaller energy 
is involved (Eq. Ill|) . 



VI. THE ROLE OF PRESSURE 

Experiments have demonstrated that hemifusion 
may be caused b y sufficiently large normal pressure 
llHelm et all I198S ) or by negative pressure i n the water 
laver ijMacDonaldL Il985l lYang and HuaneL |2002). Wc 
shall now determine the conditions under which pressure 
induced in region B can in and of itself (i.e. with no 
surface tension effects as induced by the added polymer) 
cause hemifusion expansion by forcing water to flow out 
of the contact zone. We do this by using the simplifying 
assumption that the water in region B is under a constant 
pressure p w = p n + p , where p n is the normal pressure 
on the bilayers and p is the osmotic pressure induced by 
the solute in the bulk. The finite thickness of the wa- 
ter layer in region B (whose thickness is on the order of 
a nanometer) is always maintained because of hydration 
forces: the water molecules are organized around the po- 
lar head groups of the lipids in order to partially cancel 
their electric dipole; removing the water layer would in- 
crease the free energy because of the energetic cost of 



3 soluble N-ethylmaleimide-sensitive factor-attachment protein re- 
ceptors 
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these electric dipoles whose normal components, in gen- 
eral, point to the same direction due to the hydrophobic 
nature of the lipid layer. Thus the water flow out of re- 
gion B and into region A is possible only by the expansion 
of region C. 

The energy (per unit area) difference associated with 
a pressure difference of p w is p w d, where d is the distance 
between the two proximal monolayers. This should be 
compared with the energy difference Sg associated with 
the free energy gradient i n the monolayer. In the experi- 
ment of lKuhl et al.ll)l996|) that yield p w d w 0.08 erg/cm 2 , 
which is of the same order of 8g. Nevertheless, we show 
below that the external normal pressure has only a mi- 
nor effect on the pressure in the monolayer and on its 
density. We will thus show tha t under the experimental 
conditions of IKuhl et alJ l|l996j) . the external pressure is 
insufficient to cause h emifusion expansi on. 

In the experiment o flKuhl et alJ l)l996j) the applied nor- 
mal pressure is p n = 0.3 atm and the osmotic pressure is 
p ~ 0.1 atm, so the total pressure between the bilayers is 
p w Ri 0.4 atm. We now estimate the contribution of this 
pressure to the lipid density variation in the experiment. 
For a fluid membrane, the relation between the tension 
p - the two dimensional pressure in the membrane - to 
the three dimensional pressure p w , is p w = p/h, where h 
is the thickness of the monolayer. In order to induce the 
observed density variation pj, = 0.05 the tension needed 
is \p\ = 2.5 erg/cm 2 . For h = 5 nm the pressure required 
to induce such tension is 5 atm - much larger than the ac- 
tual pressure in the experiment. Thus, the contribution 
of the normal and the osmotic pressures to the density 
variation is negligible compared with the surface tension 
effects due to the PEG-lipid interactions that result in 
densification of the lipids. This result underscores the 
point made in Sec. ]V\ changes in the pressure are much 
less effective than surface energy variation for the initia- 
tion of stalk expansion. 

We now estimate the pressure p w needed to initiate 
hemifusion, without a lipid density gradient (that is, with 
Pb = 0) . The radial force per unit length on the boundary 
at R due to the external normal pressure is 



-1 d{ Pw V) 
2ttR dR 



(13) 



From Eq. ^ the condition for spontaneous fusion is 
Pwd > ft/Ro- For the values given above, we require 
Pw > 10 7 dyne/cm 2 = 10 atm. Experimental results in 
different conditions are within that range. The pressure 
needed for the hemifusion of bilayers directly supported 
on mica (with no added polymer or other mechanisms 
that give ri s e to l ipid density gradients) w as found by 
iHelm et all l|l989l) to be p w ~ 40 atm. IWong et all 
(1999) used a surface forces apparatus to apply pressure 
on DMPC bilayers supported on polymer layers. The 
polymer layer allowed the bilayers some lateral conforma- 
tional freedom, thus permitting more freedom for the ad- 
justment of stalk shape and size Kozlovskv and Kozlovl 
(20021 iMarkin and Albanesl l2002|) . In that case, where 



the stalk geometry could easily adjust, the cost for form- 
ing the stalk was reduced and hemifusion was observed 
at a mu ch lower pressure of p w = 2 atm. In the experi- 
ment of lKuhl et al the pressure p w « 0.4 atm is 
too low to be the driving force for hemifusion. 

Pressure in itself is not enou gh to cause hem i fusion 
but it is sometimes necessary 4 . iLeckband et alJ l)l993f) 
showed that the amount of pressure needed for hemifu- 
sion is directly related to the lipid density near the con- 
tact area. In that experiment, two bilayers were brought 
into contact using a surface forces apparatus. When 
Ca ++ ions were introduced, there was a phase separa- 
tion in the bilayers. The density of lipids in the bilayer 
regions that were brought into contact was characterized 
by the hydrophobic adhesion energy. When thinner re- 
gions were brought together (characterized by adhesion 
energy of = 3.8 erg/cm 2 ) they either hemifused spon- 
taneously, or required only a small amount of pressure 
(Pn < 1 atm) to induce hemifusion. For denser bilay- 
ers (-E a d = 0.15 erg/cm 2 ) a pressure of p n = 4 atm was 
re quired for hemifu si on. 

lYang and Huan j l)2002|) induced negative osmotic 
pressure on the water layer between the bilayers by lower- 
ing the relative humidity of the environment of DPhPC D . 
At 80% humidity the lipids were at the lamellar phase. 
As the relative humidity was decreased the water were 
expelled from between the bilayers by the osmotic pres- 
sure and the lamella were connected by stalks, directly 
observed by x-ray diffraction. In this experiment the de- 
hydration was due to negative pressure of the water layer 
induced by the reduced relative humidity, and not by nor- 
mal pressure, but the physical effect of the two is similar. 



VII. SUMMARY 

In this paper we used a model based on lipid density 
gradients induced by surface energy variation that occur 
far from the hemifusion zone, to predict the the condi- 
tions for the initiation of hemifusion by stalk expansion 
and the dynamics of mesoscopic hem ifusion. Our theor y 
was motivated by the experiments of lKuhl et alJ l)l996f) . 
However, the quantitative scheme presented here can be 
generalized to any system of two lipid bilayers initially 
connected by a stalk, where a perturbation in region A, 
mesoscopically far from the stalk, causes tension in the 
membrane in that region. For example one can apply 
our results to tension induced b y the electrostatic inter - 
actions caused by calcium ions llLeckband et all Il993l) 
tensio n induced by laser tweezers l|Bar-Ziv^mdMosea. 
119941: iMoroz et all Il99fift . or the effective tension in- 
duced by the attraction of oppositely charged bilayers 



4 In the experiments it is difficult to distinguish between applied 
pressure and time in contact effects (T. Kuhl, private communi- 
cation). 
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tPantazatos and MacDonaldi ri999). 

We have compared the effect of the friction of the 
two monolayers, the the water viscosity and the intra- 
monolayer viscosity on the two dimensional lipid motion 
and showed that the friction dominates. Thus, the lipid 
dynamics depend on the friction and not on hydrody- 
namics. This means that the spacing between the two 
layers is irrelevant for the lipid dynamics. 

Experiments similar to those of iKuhl et alJ l)l996j) can 
test the predictions of the model for the time scales as 
functions of the lipid density and friction as well as the 
value of the driving force due to the tension induced in 
region A. One could vary each of the parameters pb (the 
relative change in lipid density) , a (related to the induced 
tension) and b (the interlayer friction) independently, and 
measure the 'hemifusion radius' R(t), the final radius Rb 
and the time to complete the process At as functions of 
these parameters. 

In particular, the friction b can be varied independently 
of a by changing the composition of the distal bilayers 
while maintaining the same composition of the proximal 
bilayers. The friction can be varied by changing the in- 
teractions between the chains that is responsible for most 
of the fri ction, via chain length changes or temperature 
changes IjYoshizawa et all [1993). 

Once an empirical, temporal profile for the hemifusion 
expansion, R(t), is measured for systems with known pa- 
rameters, one can use the same experiment to estimate 
the effective diffusion constant for the lipid flow, a/b, 
for different lipid bilayers. One can easily vary the lipid 
density at the boundary, <Jb by changing the polymer (or 
calcium ions) concentration since the density Ob is deter- 
mined by the equality of the chemical potentials of the 
lipids exposed to the polymer and those exposed only to 
the water. 

The static part of our theory deals with the initial 
conditions required for stalk expansion. We have eval- 
uated the necessary density variation pb = —f t /aRo and 
demonstrated that the related surface energy \ctpb 2 is 
much smaller then the surface tension apb- This result 
is not surprising, since it is a general result of a first 
order expansion around an energetic minimum. Still, 
it does give a new insight regarding biological fusion 
mechanisms. It suggest mechanisms working through the 
change of the surface energy Sg are much more effective 
than mechanisms that exert force or normal pressure on 
the stalk. 

The predicted dependence of stalk expansion on the 
lipid density can be tested by measuring the critical den- 
sity pb at which stalk expansion occurs. The results may 
serve to learn more about the stalk structure and ener- 
getics. 

We expect that near the end of the process of hemi- 
fusion expansion, when R(t) « Rb experimental re- 
sults may differ from our predictions, since the density 
profile of the polymer (or calcium ions in the case of 
iLeckband etafl Jl993)) may vary in a gradual manner 
around Rb] in our theory we assumed a sharp ('step 



function') decrease of the polymer density at Rb- We 
also expect a deviation from our theory when the ra- 
dius R(t) of the hemifusion region is close to its initial, 
molecular stalk radius Rq, due to microscopic details of 
the lipid structure in the stalk. 

We distinguish between hemifusion induced by surface 
tension gradients, which we consider in our model, and 
hemifusion induced by pressure. Hemifusio n may be in- 
duced by normal p ressur e on the bilayers ijHelm et all 
Il989t IWong et all Il999j) or by dehydration which in- 
duces negative pressure in the water layer between them 
l(Yang and HuaneL 120021) . We showed that this pathway 
to hemifusion requires much more energy (per unit area) 
than fusion that is induced by surface tension gradients. 

We have sh own that the indu ced pressure p w in the 
experiment of IKuhl et alJ (^996) cannot be the primary 
direct cause of hemifusion. Still, pressure does play an 
important role in stalk formation. It may also effect stalk 
expansion through its effect on the lipid tilt energy f t and 
on the initial stalk radius Rq. 



APPENDIX A: Local lipid dynamics 

We bring here the full calculation of the local lipid 
dynamics. Note that though in our final result we leave 
only the terms linear in p, one may also calculate in the 
same framework the non linear terms in the case p is not 
small. 

The force balance equation is 

-Vp-bv = 0, (Al) 
and the continuity equation is 

. (A2) 
Writing the energy per lipid as p(o) the surface tension 



p 



d(Np(a)) 



dA 



(A3) 



TV 



where A is the a macroscopic area and N = a A is the 
number of lipids in this area. 

From Eos. lATllAH andlAll we have 



b^- = V(aVp) = (2ctV (ct) + <r 3 p"(<j)) V 2 a 

at 



V( 

+ 5aV» + a W)) (V<t) 2 . (A4) 



To first order in the density variation p, Eq. IA4I has the 
form 



(A5) 



where a = ctq A t "( <T o) 
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APPENDIX B: Global lipid dynamics 



In Sec. IIVI we consider the boundary conditions for 
the lipid density. In order to fully predict the dynamics 
of hemifusion expansion, we also need to determine the 
flow at the boundaries. For this we use the integral form 
of the continuity equation: 



d [ Rb [ _ - 

— I 2irrdr air) = — d> av ■ dl . 
dt J R J Rb 



(Bl) 



The left hand side of Eq. |BT] describes the rate of change 
of the lipid number in region B while the right hand side 
gives the flow of lipids through the boundary Rb- We 
assume cylindrical symmetry, so v = v r f. From Eqs.EU 
and lA3l we obtain 

Vr (r) = -l(2 f j,'(a)+*n"(*))a^-. (B2) 

We now use Eq. IA5l to calculate the left hand side of 
Eq. |B1] 



dt 



2-Krdr er(r) 



„ a I _ da 

2 'i ( R °i» 



or 



If we take only terms linear in p, Eq. IB II gives: 



/ Rb dp 




dp 




I 2R dr 


Rb 


dr 





-2vRa{R)^. (B3) 



(B4) 



APPENDIX C: Adiabatic solution 

Equations IB4I and IA5I along with the boundary con- 
ditions completely determine the time evolution of the 
monolayers to first order in p. From these equations we 
can calculate R(t), and predict the temporal profile of 
hemifusion expansion. We write these two equations us- 
ing dimcnsionless variables and scale the spatial variables 
so that they are of order of unity, in order to get an esti- 
mate of the time scales. The natural spatial scale is the 
final size of the hemifusion region, Rb- We thus define: 
x = r/Rs, R = R/Rb-, P = p/pb as well as two time 
variables: a "fast" time r = at/bRs 2 at which the local 
lipid flow occurs, and a "slow" (since pi, is small) time 
f = \pb\r which is the scale over which the hemifusion 
region expands. Eas. IB4I and I A5I become 



dR 

df 



dp 
dr~ = 

= d -l 

dx 



op 

dx 2 



x=R 



_ 

x dx ' 
2R dx 



(CI) 



(C2) 



Since all the variables that appear on the right hand 
side of Eqs. ICll and lC2l are of order unity, the units of r 



and f suggest the time scales of the processes described 
by the equations. For \pj,\ <C 1 we have r ^> f , which 
implies that we can use an adiabatic approximation: the 
local lipid flow occurs quickly so that the lipid density 
is instantaneously given by the asymptotic equilibrium 
solution of Eq. ICll for r — > oo. We then use this solution 
to determine the slower time evolution of the hemifusion 
radius R from Eq. IC2I 

At asymptotically long times, both sides of Eq. IC1I 
vanish. The adiabatic density profile reached is 



p(x) = 1 - 



log(x) 



(C3) 



(C4) 



log(fl) 

Plugging this solution into Eq. IC2l we obtain 

dR _ -1 
df ~ 2R\og(R) ' 

The solution of this equation is implicitly given by 

2f = R 2 {1 - 21og J R) . (C5) 
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